UAB-FT-725 
LPT-ORSAY/12-81 



Implications from clean observables for the binned 
analysis of B — > K* ji^ at large recoil 

Sebastien Descotes-Genon 

Laboratoire de Physique Theorique, CNRS/Univ. Paris-Sud 11 (UMR 8627) 

91405 Orsay Cedex, France 

JoAQUiM Matias, Marc Ramon, Javier Virto 

Universitat Autdnoma de Barcelona, 08193 Bellaterra, Barcelona, Spain 

Abstract 

We perform an extensive analysis of g^-dependent B — )• )■ K'it)£'^£~ angular 
observables at large recoil, aiming at bridging the gap between current theoretical 
analyses and the actual experimental measurements. We focus on the most appro- 
priate set of observables to measure and on the role of the g'^-binning. We highlight 
the importance of the observables Pi exhibiting a limited sensitivity to soft form 
factors for the search for New Physics contributions. We compute predictions for 
these binned observables in the Standard Model, and we compare them with their 
experimental determination extracted from recent LHCb data. Analysing 6 — )• s 
and b — )• st^t~ transitions within four different New Physics scenarios, we identify 
several New Physics benchmark points which can be discriminated through the mea- 
surement of Pi observables with a fine (/^-binning. We emphasise the importance 
(and risks) of using observables with (un) suppressed dependence on soft form factors 
for the search of New Physics, which we illustrate by the different size of hadronic 
uncertainties attached to two related observables (Pi and S's). We illustrate how 
the g^-dependent angular observables measured in several bins can help to unravel 
New Physics contributions to B ^ K*{-^ KTT)i'^i~, and show the extraordinary 
constraining power that the clean observables will have in the near future. We pro- 
vide semi-numerical expressions for these observables as functions of the relevant 
Wilson coefficients at the low scale. 



1 Introduction 



The set of rare 5-meson decays mediated by the 6 — )■ s transition has been thoroughly 
studied for many years both from the theoretical and experimental sides. Recently, this 
set has been complemented by increasingly precise experimental measurements of q^- 
dependent angular observables in the decay B — )• > K7r)£+£~ [H El El HI |5]. This 
has triggered extensive theoretical work studying the constraining power of radiative and 
semileptonic B decays on New Physics in the framework of the weak effective Hamilto- 
nian [El El [HI [HI [TOl [HI [12]. These constraints apply mostly to the Wilson coefficients 
C7 , Cg and C\q related to the magnetic and semileptonic operators Oj, Og, Oiq and 
the corresponding chirality-flipped operators (which are highly suppressed in the Stan- 
dard Model). In addition, the tight experimental bounds set in the last few months on 
the differential decay rate of Bg — )■ [111 [HI [151 [Ej, pushed it close to its small 

Standard Model (SM) prediction [TTl [T8] - the theoretical and experimental values are 
actually brought even closer by the Bg mixing correction ©(AF^) to Bs branching ratios 

measured at LHCb, discussed in Refs. [HI [20] and applied to Bs — )■ in Ref. pT] . 

(') (') (') 

This puts strong constraints on C)q , as well as on the coefficients Cy and Cp of scalar 

and pseudoscalar operators, specially when combined with B — )■ Ki~^i~ data |22l [23] . 

Small experimental errors and a good control over hadronic uncertainties on the theory 

side are the key ingredients for these constraints to be efficient. At present, the branching 

ratio of the inclusive radiative decay B — t- XsJ, and the CP asymmetry of -B — ?■ K*'~f 

constitute the strongest constraints on the C7-C7 plane. However, the complementarity 

of constraints among different observables can be exploited to reduce considerably the 

parameter space. In this respect, the inclusion of the isospin asymmetry of 5 — )■ K*'~f, 

together with the forward-backward asymmetry Aps and the longitudinal polarization 

fraction Fl in B ^ K*i~^i~ integrated over the dilepton invariant mass between 1 and 

6 GeV^ already impose additional nontrivial constraints on the C7-C7 plane, as well as on 

cP andcS PE]. 

This complementarity can be exploited further by considering the g^-dependence of 
angular B — > K*i'^i~ observables. Indeed, increasingly precise measurements of these 
observables integrated in smaller bins are being presented, due in part to the important 
statistics obtained from the large data sets collected at the LHCb experiment. This in turn 
allows for more complete angular analyses providing more observables |1] . The prospects 
for the near future are very good, aiming towards a complete angular analysis with a fine 
g^-binning. 

At this point, it becomes crucial to handle theoretical uncertainties as accurately as 
possible. The lack of huge deviations in B physics up to now forces us to be precise and 
conservative in our theoretical predictions. Before claiming any discrepancy, indication 
or discovery of New Physics, one must be sure that a description in terms of observables 
with little dependence on the specific choice of hadronic parameters has been used. While 
the uncertainties related to the form factors constitute a dominant part of the theoretical 
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error, there is a wide spread of quoted uncertainties for B — > K* form factors in the recent 
hterature, ranging from a ~ 10% to a ~ 40% error for the same form factor [211 [25] . 
The size of this error does not depend only on the particular theoretical method used 
to compute the form factor, but also on the delicate estimation of errors associated to 
the assumptions built in each procedure. For example, the values Aq{0) = 0.33 ± 0.03 
and ^(0) = 0.31 ± 0.04 given in Ref. [2lj should be compared to the values Ao{0) = 
0.29 ±0.10 and V^(0) = 0.36 ±0.17 as quoted in Ref. |25]. Even central values have shifted 
significantly, see for instance the value V{0) = 0.41 ±0.05 from Ref. [26] before its update 
of Ref. [21] (also consistent with Ref. [2?]). Without attempting to discuss the related 
conceptual issues in any further depth, it is clear that the impact of such discrepancies 
between different groups concerning hadronic uncertainties can be reduced greatly if one 
selects quantities that show a suppressed dependence on the soft form factors, such as the 
observables or A^j, defined in Refs. [23 1211 [30] . These observables can be considered as 
being theoretically clean in the kinematic range of interest here. 

The construction of theoretically clean observables in B K*i~^i~' has been the 
subject of theoretical work for some time [211 [291 [311 [321 [301 [33]. Based on the symmetries 
of the B — )• K*i'^i~ angular distribution discussed in Ref. [31], a complete characterization 
of the full distribution in terms of a minimal basis of clean observables has recently been 
proposed [30]. The conclusion is that, a complete description of the differential decay 
rate in the limit of massless leptons and in the absence of scalar contributions can be 
achieved through a set of 6 clean observables -Pi,2,3,4,5,6 complemented by two observables 
with a significant sensitivity to form factors, e.g., the differential decay rate and the 
forward-backward asymmetry ApB (or equivalently F^) 0. In a short term, these clean 
observables (Pj) should play a leading role in detecting deviations from the SM in a safe 
way, relegating less clean observables such as (Afb), (-^l) or {S3) to a secondary role of 
useful cross-checks. 

(2) 

Measurements of the transverse asymmetry Aj, have been already provided by the 
CDF collaboration [2] . However, there has recently been some reluctance from the exper- 
imental side to extract such clean observables because their theoretical predictions have 
been mostly presented as functions of g^, while the experimental results are obtained 
integrated in g^-bins. This issue becomes relevant when a coefficient in the angular dis- 
tribution is expressed as a product of various observables. For example, in Ref. [4], the 
authors prefer to fit for and 5*3 instead of F^ and Aj.' (where 5*3 ~ (1 — Fl)Aj.'), argu- 
ing that a rapid variation of both F^ and A^"* with could result in a biased estimate of 
A'^^ when averaging over large g^-bins. This is a perfectly correct statement if one wishes 
to compare the experimental measurements with (differential) g^-dependent observables 
such as A'^\q'^). Actually, in the specific case of A^^ (for the SM case as well as for 

^Mass and scalar effects can be taken into account by including 4 extra clean observables (see Ref. [30]). 
The explicit form of the symmetries in the presence of mass terms and scalar operators can be found in 
the same paper. 
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particular NP values of the Wilson coefficients involved) , we observed a very good agree- 
ment between the observable integrated naively over the low-g^ bins (i.e., J^.^^dq'^Aip^) 
and its value derived from the binned observables as they are being measured (denoted 

(2) 

(y4^ )bin and defined in the next section), showing that the bias is small in this partic- 
ular case. Though encouraging, this remark will not prevent us from trying to simplify 
the comparison between theory and experiment by providing theoretical predictions for 
the exact integrated quantities as measured in experiments. These kind of integrated 
observables have been already discussed in the context of other transverse asymmetries 
in Refs. [331 [11]. 

In the present paper we address precisely how to analyse efficiently the LHCb mea- 
surements on 5 — )■ K*i^i~ at large recoil by choosing a set of clean observables integrated 
over g^-bins. We begin by building a set of integrated observables that correspond in the 
limit of small binning to the observables in Ref. [30], with the aim of making contact with 
experimentally measurable quantities. In Sections [2] and |3l we present Standard Model 
predictions for these integrated observables, providing the results for different choices of 

the g^-binning. In Section H] we perform a model-independent analysis setting constraints 

(') (^) (') 

on the Wilson coefficients Q , Cg and CIq, using data from radiative B decays and in- 
cluding the forward-backward asymmetry and the longitudinal polarization fraction of 
B — )■ K*i~^i~, both integrated in the full low-g^ bin [1, 6] GeV^. This updates the analy- 
sis of Refs. [6l [7] (and related analyses) with several improvements. We use these results 
to identify a number of New Physics "benchmark points" that are allowed at the 95.5% 
confidence level by all the constraints considered. In Section [5] we discuss the potential 
complementarity of g^-dependent observables in i? — )■ K*i'^i~ by analyzing our set of 
clean observables within the NP scenarios specified by the benchmark points. This analy- 
sis indicates the scenarios that are more likely to be affected predominantly by the binned 
observables in _B — )■ K*i^i~. 

In Section [6] we open up a parenthesis to discuss the impact of hadronic uncertainties 
on the different observables, and the resulting model-independent constraints that follow 
from different choices of observables. We demonstrate the advantages of using a complete 
set of clean observables as the one introduced in Section [2] and Ref. [30] . 

We then address the model-independent constraints from g^-dependent observables in 
Section [71 We extract the experimental values for the three clean observables (Pi, P2 and 
P3) that can be related to the measurements provided in Ref. Our determinations are 
affected by uncertainties that could be considerably improved, since we lack the experi- 
mental information concerning the correlation between the LHCb measurements, which 
is essential to assess uncertainties in a proper and accurate way. The constraints from 
g^-dependent observables are studied in Section [H] First we consider the constraints from 
the measured ApB and in the two g^-bins [2,4.3], [4.3,8.68] GeV^, then we turn to 
the constraints imposed by the clean observables Pi, 2,3 in the same bins. In Section [9l we 
compare briefly our study with other similar works. 
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After concluding in Section [TOl we include in Appendix [A] the expression of the coef- 
ficient Js in term of observables and in Appendix |B] a collection of the relevant formulas 
used to derive the New Physics constraints. In particular we provide the explicit ex- 
pressions of the clean integrated observables used throughout the paper, for the different 
choices of g^-binning. In Appendix [C] we describe the statistical approach used in the fits. 

As a summary of the most important findings, we anticipate the following conclusions 
of this work: 

• Three (Pi, 2,3) out of the six clean observables describing the massless distribution can 
already be extracted from current measurements, as shown in Section [71 While P3 
sets no relevant constraints yet. Pi and P2 are complementary to other radiative and 
semileptonic observables. Combining the measurements of P2 in two different bins 
leads to a mild tension with respect to the SM, marginally compatible at 95.5%C.L. 
and pointing towards a negative contribution to the Wilson coefficient C7. 

• The explicit form of the coefficients of the massless angular distribution in terms of 
the basis of observables is given in Eq. (HM . The expressions turn out to be very 
simple and exhibit two important features: 

— A more natural basis devised to extract information from the distribution in 
a clean way emerges in the massless case, with a slight redefinition of the 
observables ^4,5,6 — ^ P456' which are also clean and defined in Eqs. (IT5l)-(fT71) 
(in the SM ^4,5,6 ~ P4 5 e ^ ^^^y good approximation). 

— The clean observables in the natural basis (Pi, 2,3 — ^ -^4,5,6) "^^^ be related to 
form-factor sensitive observables Si through the following simple rule: 

^3fi^ p '^4,5,7 ^ p/ 

where the exact relationships are given in Eqs. ([5]) and ( ITSll . 

• The "fiipped-sign solution" for C7 is in general disfavoured by present data at the 
95.5% confidence level depending on the NP scenario considered. The isospin asym- 
metry in P — )■ K*'~f plays an important role (independent of Wilson coefficients 
other than Q ) , as well as the forward-backward asymmetry in P — )■ K*i£. This 
confirms the result of Refs. [HI [7]. 

• We show explicitly the strong impact of the different computations available in the 
literature for the soft form factors on the theoretical uncertainties for observables like 
ApB, Fi and 5*3, and the robustness of the clean observables Pj. While the impact 
on the theoretical error in Fl is evident, the problem for observables like is more 
subtle. In the case of 5*3 the theoretical uncertainty in the SM is protected by its tiny 
central value, but away from the SM point the impact can be substantial, preventing 
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this observable from discriminating NP scenarios. None of these problems affect the 
clean observables Pi or Alp. 



2 Integrated observables in g^^-bins 

The differential decay rate of the process — )■ K*{-^ Kti)['^1'~ can be written as: 
d^T 9 r 

— — -— — — = — — Jis sin^ Ok + -he cos^ Ok + ( sin^ Ok + J2c cos^ 9k) cos 26'; 

ag^ ttcos Uk acos Ui dcp SZrr 

+ J3 sin^ 9k sin^ 9i cos 20 + J4 sin 26^7^ sin 29i cos (p + J5 sin 26*7^ sin 6'/ cos (p 

+ {Jqs sin^ 6'/^ + Jqc cos^ 6'x) cos 9i + J7 sin 26'i^ sin 9i sin + Js sin 29k sin 26*/ sin 



+ Jg sin^ 6'x sin^ 6'; sin 20 



(2) 



where the kinematical variables 0, ^i^, are defined as in Refs. [331 EH EO] . The decay 
rate f of the CP-conjugated process Bd — > K*{-^ KiT)i~^i~ is obtained from Eq. (E]) by 
replacing Ji, 2,3,4,7 — ^ -^1,2,3,4,7 and J5,6,8,9 —-^5,6,8,9, where J is equal to J with all weak 
phases conjugated. This corresponds to the same definition of 9£ for both B and B (see 
for example [SUElj). In this paper we assume that all the observables are CP- averaged, 
and so are always functions of Ji + Ji. Therefore, Ji ^ Ji + Ji and F — )> F + f should 
be understood in all the formulas below, and in particular all the observables 0{J) are 
assumed to be 0( J + J). 

In order to cope with limited statistics, one can write down integrated distributions, 

such as the uniangular distributions, which depend on a subset of coefficients Ji. This is 

(2) 

the way observables such as Fl, ApB or Aj, have been measured traditionally. A more 
recent approach deals with "folded" distributions, with the double advantage of increasing 
the statistics and focusing on a restricted set of angular coefficients. For example, in 
Ref. ji], the identification of o + tt has been used to produce a "folded" angle 
G [0,7r] in terms of which a (folded) differential rate dr{(f)) = (iF(0) + (iF(0 — tt) 
becomes 
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dq"^ dcos 9k dcos 9i 



IGtt 



Jic cos^ 9k + Jisi^ — cos^ 9k) + J2C cos^ 6'x(2 cos^ 9i — 1 



) 



+J2,(1 - cos'^9K)i2cos'^9e - 1) + ^3(1 - cos^ 9k){1 - cos^^^) cos: 



-|- J6s(l — cos^ 9k) cos9i + 79(1 — cos^ 6'x)(l — cos^ 9i) sin 20 



(3) 



Neglecting scalar and mass effects, this distribution can be written as a function of the 
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observables in Ref. [30] as follows: 



d^T 9 



dq'^ dcos 6k dcos 6i d(j) ISvr 



Fl cos^ + ^Ft(1 - cos^ Ok) - Fl cos^ ^^^(2 cos^ e^-l) 



+-Ft{i - cos^ eK){2 cos' et-i) + ]:PiFt{i - cos' eK){i — COS 6'^) COS 20 



+2P2i^T(l - cos' ^i^) COS 6*^ — P^Ft{^ — cos' 6'i<')(l — cos' 6i) sin! 



— (4) 



where Pi, P2 and P3 are theoretically clean observables ^30j that in terms of form factor 
dependent observables are given as 

PiPt = 2S'3, P2Pr = ^6s/2, PgPr = -^9 , (5) 

or alternatively, 

P^Ft = 2Afb/3 , P3PT = -Am . (6) 

The quantity Ft is defined as Ft = ^ — F^. The observable Pi is better know by its 
original name, Aj, [28] . 

Experimentally, one can fit separately each of the following five independent coeffi- 
cients that appear in the partial distributions: 



(7) 

For each g' one can then in principle extract the theoretically clean observables Pi = 
Aj.\, P2 and P3, as well as the transverse polarization fraction and the differential decay 
rate. However, in practice the g'-dependence is discretised in a number of bins, and the 
coefficients Cj(g') thus extracted are quantities integrated over particular g'-intervals: 

_ /bin^?'c,(g') 

\'^i/bin 



One could hope to be able to extract also corresponding integrated theoretically clean 
observables such as 

/ ^ (2) \ (naive) 

but due to the experimental procedure used, such a determination is achievable asymptot- 
ically only, as the bin size goes to zero. The actual theoretically clean quantities that can 
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be extracted from experiment and on which we will focus from now on, must be composed 
of the integrated quantities in Eq. (|8]): 



\A/bin = /bin - 7 - 1 T , UUj 

Jbin^^C4(g2) (C4)bin 

/p^ _ /bin^g'c2(g^) _ (C2)bin , . 

Jbin^^C4(g2) (C4)bin 

bin / 1 (->\ 

Other observables are accessible to the current LHCb data set by means of similar partial 
angular analyses. Three observables related to P4, P5 and Pg (see Ref. [30]) could be 
extracted by the LHCb collaboration in the near future. This means that, without actu- 
ally performing a full angular analysis, the LHCb collaboration could be able to provide 
measurements of the complete set of 8 observables that describe the full distribution in 
the massless approximation (six of them being theoretically clean [30]). Therefore, we will 
work here under the assumption that integrated versions (exactly as in Eq. ([8])) of the 
observables Co,i, 2,3,4 and ^4,5, 7 are available experimentally. Any measurable observable 
must be a combination of the (cj)bin in Eq- (E]), and of the observables: 

/J ^ _ibin^^jW^ 

('^4,5,7)bin - J-^, . (13) 

The coefficients of the angular distribution can be written in terms of the basis of 
observables (see Ref. [30]), and in terms of the coefficients q of Eq. ([7j) as follows: 

3 dV 3 1 rfr 1 



dV dV 

^L-r^ — Co — C4 , J2c — —-TL-r^ - 
dq'^ dq"^ 

Vp -V 7 -2PP 



J5 = P^^F^L ^ = ^5^04(00-04) , 

J7 = -Pl^^F^L ^ = -^6^04(00 - C4) , (14) 
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where the primed observables are defined as: 



Pi^P,V^^i= (15) 

H^p,v^n\ = -jd^ (16) 

■^V ~'J2c'J2s 
J7 



P',^P,^l-P, = --^J== (17) 

The case of the coefficient Jg is discussed separately in detail in Appendix |X1 The ex- 
pressions (fT5|) -( IT7|) in terms of Jj are exact even in the presence of scalar contributions 
and mass terms. These observables -P^ 5 g are clean and coincide to a good approximation 
with -P4,5,6 in the SM (due to the fact that Pi in the SM). The whole analysis can 
be performed directly in terms of the observables -P4,5,6 i however, from the experimental 
point of view, fitting the primed observables is simpler and more efficient. 
These observables can be related to the observables 6*4^5 j of Ref. [2^ : 

P'^ = 2 , = , = (18) 

vFtFl vFtFl a/ FtFl 

There is therefore no particular advantage for the experimental extraction of the observ- 
ables Si instead of the P/, while from the theory point of view the P/ are under better 
control and suffer from smaller uncertainties. 

We want now to construct the theoretically clean integrated observables that corre- 
spond to those in Ref. [20] (or variations thereof). For Pi, 2,3 the answer is precisely (Pi)bin5 
(-P2)bin and (P3)bin dcffucd in Eqs. fll0m2p . In analogy with Eqs. (11011121) . integrated ver- 
sions of the observables P^ 5 g can be deffned: 



/p/\ Jbin"-'^ -^syj ; yj^ihin 

(Ps/bin = , = ; = , (20) 





- c,{q^)) ^ 


\J (c4)bin((Co)bin ~ (C4)bin) 






(^5)bin 


\l /bin dq'ci^q') /bi, dq^{co{q^) 


- c,{q')) ^ 


\J (c4)bin((Co)bin " (C4)bin) 


-Lndq'Jriq') 




— (-^7)bin 



(^^)bin=^= ^ ^— = ^ ^ (21) 

y Ihindq^C^iq^) Iundq^iCo{q'^) - C^iq^)) y (c4)bin((Co)bin - (C4)bin) 

Finally, integrated versions of the longitudinal polarization fraction and the forward- 
backward asymmetry Ap^ can be deffned in terms of the coefficients Cj in the following 
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g2 (GeV2) 


[1,2] 


[2, 4.3] 


[4.3, 6] 


[1,6] 


(^i> 


n nns+O-009+0.051 

U.UUO_Q QQ5_Q 


n nc;i +0.007+0.046 
^•UO±_Q oii_o.051 


n nQQ+O-OOl+0.049 

u.uyy_o 001-0.053 


n nc:.t:+0.007+0.048 
U.UOd_o 009-0.053 


(^2) 


n 009+0.023+0.013 
"-'•"-'^^-0.013-0.01 


PI 909+0.032+0. 
'-'•^"'^-0.105-0.032 


r, 9 cr +0.048+0.025 
'-'•^^-0.081-0.045 


n ns/1 +0-035+0.008 

U.UO'i_o 092-0.034 


(^3) 


U.UU^_Q QQ3_Q Q2g 


U.UU'±_Q 004-0.025 


n nn'?+0-00i+0.023 

U.UUO_o 002-0.024 


n nn'i+0-+o.023 

U.UUO_o 003-0.025 


(Pi) 


n 1 70+0.022+0.026 

'-'•J- '^-0.008-0.025 


n Kcfi+0. 068+0.008 
U.OOD_o. 031-0. 


n qQq+0.024+0. 
u.yoy_o.oi8-o.oo3 


n c:oQ+0. 061+0.01 
U.00y_o. 027-0.001 


{Pi) 


U.OUU_Q Qg7_Q Q2g 


n o/i 9+0.056+0.022 
U.O4:^_0 i44_o.07i 


n 77c;+0.047+0.08 
''-'-0.072-0.106 


n q/iQ+0.052+0.032 
U.O^y_0 125-0.075 


(Pi) 


n 1 09+0.008+0.005 


n ini +0-011+0.02 

U.lUl_o. 047-0.029 


n 089+0-012+0.045 
U.UOZ_o o4_o.051 


n nQ'?+0-Oii+o.026 

U.UJa_o.044-0.034 



Table 1: SM predictions for the clean observables (Pi). 



way: 

_ /bin^g'(go(9') - C4(g')) _ (Co)bin " (C4)bin . . 

^ Ln^^^^o(g^) " (Co)bin • ^''^ 

In the following sections we will study these integrated observables in detail. 



3 SM predictions for integrated observables 

We can provide SM predictions for the set of integrated observables (Pj) as well as (Afb) 
and (Fl). In Tables [T] and [2] we show the predictions in the g^-bins [1,2], [2,4.3], and 
[4.3,6] (GeV^) - following the binning used by the experimental collaborations up to now 
(except for the first bin) - as well as the predictions for the integrated low-g^ observables, 
in the region [1,6] GeV^. The first error accounts for all parametric uncertainties, while 
the second error corresponds to an estimate of A/rrib corrections, as described below. In 
Figs, m [2] and [3] we show the corresponding SM predictions for the observables in the case 
of one and three bins (corresponding to the predictions in Tables [T] and [2]) , as well as for 
five bins with a width of 1 GeV^. 

The SM predictions are obtained as follows. The observables integrated over each bin 
are defined in terms of the coefficients Ci{q^) in Eqs. ([10])- ([12]), (IIS]) -([21]), (|22]) and ([23]). 
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g2 (GeV2) 


L ' J 


[2, 4.31 


[4.3, 61 


[1,61 


(Afb) 


n 1 cr +0.041+0.007 
'-'■-'-^-0.024-0. 


n nc;Q+o.o3i+o.oo7 

U.UOO_Q oii_o. 


n n7fi+0-029+0.016 
U.um_o.oi8-0.009 


n no/i +0-027+0.01 

'-'•'-'^^-0.011-0.001 


{Fl) 


n 7/1(3+0.076+0.015 
'^"-0.039-0. 


n 0/17+0.039+0.008 
U.O^/ _0.028-0.002 


n 7Q7+0.039+0.004 
^■li" -0.033-0.004 


n QnQ+O-045+0.008 
U.OUy_Q 032-0.002 



Table 2: SM predictions for (Afb) and (Fl). 



The coefficients Cj(g^) are simple functions of transversity amplitudes (see for example 
Ref. [3DJ). The transversity amplitudes can be written in terms of Wilson coefficients and 
B — > K* form factors following Refs. [5^ Concerning the Wilson coefficients, the 
form factors, and the treatment of uncertainties, we proceed as in Refs. [501 El E] , with a 
slight revision in the treatment of form factors: 

Wilson coefRcients: The SM Wilson coefficients are evaluated at the matching scale 
/io = 2M]y, and evolved down to the hadronic scale /it = 4.8 GeV following Refs. |3Z1 ESI 
ESI HOI SI]. The running of quark masses and couplings proceeds analogously. The SM 
Wilson coefficients at the scale /if, are shown in Table El 

Form factors: There are seven B — )■ K* form factors: V(g^), ^0,1.2(9^) and Ti 2,3(1?^). 
Their determination involves the computation at = and the parameterization of the 
g^-dependence. At = 0, these form factors can be obtained from light-cone sum rules 
with B- meson distribution amplitudes (see Refs. [2S1I2S])- Concerning their dependence 
on g^, Ref. [25] provides a conservative and convenient parameterization (the prospects for 
these form factors from lattice QCD have been discussed in detail in Ref. [32] )• The soft 
form factors ^|[^^ are defined in terms of the full form factors following Ref. [SB]. The soft 
form factor ^±{0) at g^ = is obtained directly from V{0) as given in Ref. [2S] (updated in 
Ref. [25)- The form factor ^||(0) is defined as a combination of the form factors ^i(O) and 
^2(0), and in the large- recoil limit and at leading order in it is proportional to Aq{0). 
We use Aq{0) to fix ^||(g^) at g^ = to a good accuracy, and set its g^-dependence to 
reproduce its exact expression in terms of Ai(g^) and A2{q'^) using the parameterization 
of Ref. [25] . The numerical inputs used are collected in Table El 

In the present paper we take the form factors of Ref. [21] (which updates those of 
Ref. [26]) to facilitate the comparison with other recent analyses in the literature. This 
choice is different from that in Ref. [30], where the form factors of Ref. [25] were used. 
The use of the updated value for ^(0) from Ref. [21] in the present analysis shifts the 
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ruB = 5.27950 GeV [13] 


rriK* = 0.89594 GeV [43j 


'it'Sg — o.ouuo vjft: V [ I 


IIL^ — U.XUOUOOOUl vac V 


sin^ 9w = 0.2313 03] 
= 80.399 ± 0.023 GeV p] 


Mz = 91.1876 GeV [33] 


«em(Mz) = 1/128.940 i42j 


asiMz) = 0.1184 ± 0.0007 [43j 


mP°''' = 173.3 ± 1.1 GeV [S] 
fn^ \iiic) — i.zi ic u.uy Vore V [^o] 


m^^ = 4.68 ± 0.03 GeV [45j 
TTig Voce V ) — u.iui m u.uzy vjev [^<Jj 


XcKM = 0.22543 ± 0.0008 [16] 
p = 0.144 ±0.025 [16] 


Aci^^^ = 0.805 ± 0.020 [46] 
r/ = 0.342 ± 0.016 [46] 


^X^ez/) =0.1061 ±0.00017 [12] 
A2 = 0.12 GeV^ [12] 


C = 0.58 ±0.016 [42] 


Ah = 0.5 GeV [17] 
fK',\\ = 0.220 ± 0.005 GeV [21] 
\/(0) = 0.311 ±0.037 [21] 
ai,||,±(2 GeV) = 0.03 ±0.03 [21] 
AB(/i;,) = 0.51 ± 0.12 GeV [21] 


/b = 0.190 ±0.004 GeV [IE] 

/i^-,±(2 GeV) = 0.163(8) GeV [24] 
ylo(O) = 0.333 ±0.033 [21] 
a2,||,±(2 GeV) = 0.08 ±0.06 [21] 


/b^ = 0.227 ± 0.004 GeV [H] 


tb, = 1.497 ±0.015 ps [13] 




C2(/Ufe) 






CnifJ'b) 


Cel/^b) 


cf{^ib) 


cf{^^b) 


Cgif^b) 


Cio(/Wfe) 


-0.2632 


1.0111 


-0.0055 


-0.0806 


0.0004 


0.0009 


-0.2923 


-0.1663 


4.0749 


-4.3085 



Table 3: Input parameters used in the analysis and Wilson coefficients at /i;,. 
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Figure 1: Binned Standard Model predictions for the clean observables (Pi.2,3), for a single 
bin [1,6] GeV^ (left column), three bins [1,2], [2,4.3], [4.3,6] GeV^ (central column), and 
five bins of width 1 GeV^ each (right column). The red (dark gray) error bar correspond to 
the A/mft corrections, the yellow one (light gray) to the other sources of uncertainties. If 
one of the two bands is missing, it means the associated uncertainty is negligible compared 
to the dominant one. 



central values of A-p^ and in the whole low-g^ region, while only a mild effect around 
~ 6 GeV^ could be seen in some of the Pi observables. In addition, we believe that 
the impact of more recently published determinations of the form factors in Ref. [25] 
should also be considered. These form factors have much larger uncertainties than those 
of Refs. [231 and translate into large error bars in A-p-Q, Fi and other form factor de- 
pendent observables. This is a clear indication that the theoretical uncertainties attached 
to these observables (i.e., A^-q, Fi) in the literature should be considered with a healthy 
dose of skepticism, and are strongly dependent on the choice of the B — )■ K* form factors. 
Fortunately, it will be seen in Section [S] that the error bars for the clean observables P, 
are not affected by this variation of the form factor uncertainties and remain under good 
theoretical control. 
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Figure 2: Binned Standard Model predictions for the clean observables {Pi^^^o), with the 
same conventions as in Fig. [H 
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Figure 3: Binned Standard Model predictions for the observables {A-fb) and {Fl), with 
the same conventions as in Fig. [H 
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Uncertainties: We consider five main sources of uncertainties: tlie renormalization scale 
fib, tlie quark mass ratio rhc = mc/mi,, the form factors, the factor that determines the 
relative size of the hard-scattering term with respect to the form factor contribution to 
the amphtude (defined in Eq. (55) of Ref. |35]) and the A/rrib power corrections. 

We follow the usual procedure consisting in varying the renormalization scale from 
to 2fif,- For rhc we take rhc = 0.29 ± 0.02 (see Refs. |35l |19]). Concerning the 
form factors: we express all the observables as a function of Ao{q'^) and Vlq"^). The q^- 
parameterization of these form factors (that of Ref. p5]) depends on the normalisation 
parameters Aq{0), V{0) as well as and fof" (encoding the g^-dependence of the form 
factor in the so-called z-parametrisation) . We vary simultaneously ^(0) and , and in- 
dependently, Ao{0) and bf°, within the errors quoted in Tableland Ref. [25]. The two 
errors are added in quadrature. Concerning the parameter describing the relative weight 
of the hard-scattering contribution compared to the form-factor one in Ref. [35], its error 
is estimated at the level of a 25%, where its reduction with respect to Ref. [30] (where 
it was 30%) is due to the updated value of [H]. Finally, for A/ nib corrections, we 
follow the statistical procedure outlined in Ref. [31] to produce upper and lower la ranges 
consistent with a generic 10% contribution of power corrections to the amplitudes. All 
individual uncertainties are considered separately and their impact on each observable is 
monitored to produce asymmetric upper and lower errors. All upper and lower uncertain- 
ties, excluding A/rub corrections, are added separately in quadrature to produce the first 
asymmetric error bars shown in Tables [1] and [21 and the yellow bands in Figs. [H [2] and [31 
The second error bars in Tables [H and [2l correspond to the A /nib corrections. In Figs. [H 
[21 and [3l both uncertainties are added linearly to give the larger red error bands. 

As can be seen from the plots, some observables appear to be almost insensitive to 
power corrections. This merely indicates that the hadronic uncertainties are much more 
important than power corrections for these observables. 

4 Model-independent constraints without g^-binned 
observables 

In this section we revisit the model- independent constraints on the Wilson coefficients C7, 
C7, Cg ^ and C[q from well controlled observables, excluding all B — )■ K*£~^£~ observables 
except for {ApB)[ifi] and (Fl)[i,6]- This analysis follows closely the study of Ref. [6l [7] 
(see also Refs. [211 IHl M lEl [50]). The aim is to reevaluate the constraints on the Wilson 
coefficients taking into account the following updates: 

• Updated averages for {Apb)'[iq] and (-FL)[fe] including the new measurements in 
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Observable 


Experiment 


SM prediction 


BR{B X,7) 


(3.55 ±0.26) ■ 10-^ 


51 




(3.15 ± 0.23) ■ 10"^ [321 






-0.16 ±0.22 [5l] 


-0.03 ±0.01 [6j 






0.052 ±0.026 [51] 


0.041 ±0.025 [61 




BR[B Xsfi^fi )[i,6] 


(1.60 ±0.50) ■ 10-6 






(1.59 ± 0.11) • 10-6 [52 




(^Fb)[1,6] 


"-'•^'-'-0.078 ^ 


-0.024 ±0.020 t 


(i^L)[l,6] 


622+°-°^^ * 


0.809 ±0.039 t 




< 4.5 ■ 10-*^ (at 95.5% C.L.) [H] 


(3.32 ± 0.17) ■ 10-9 [12J 



Table 4: Experimental numbers and Standard Model predictions of the observables used 
in the analysis of Section |H ic indicates it is our own average of the data, f indicates it 
is our SM prediction. (See also Ref. [T7] concerning the Bg — )■ fj.~^fj.~ branching ratio.) 



Ref. [1] by the LHCb collaboration: 



/A \LHCb _ n 1 8+0-06+0.01 fr)A\ 
V^FBhi 61 - -U-iO_o.06-0.02 l^^j 



.,6] 

HC 

.,6] 

The experimental averages for these two observables are collected in Table |H 



/p \LHCb _ ^^+0.06+0.04 (r)rT\ 
\^L)h(i] - 0.66„o.o6-0.03 l^5j 



Updated theoretical predictions for (v4fb)[i,6] and (-Fl)[i,6] including subleading cor- 
rections of order iV^^V^sl/lVf^VtsI, as well as a recent update for fs [H]: 

fB = 190 ± 4 MeV (26) 



• Analysis of the constraints using a consistent (frequentist) statistical approach de- 
tailed in App. O 

These results will be used in the following sections to study the impact of the inclusion of 
binned observables, in view of (a) recent measurements at LHCb and (b) the impressive 
prospects for the near future measurements of g^-dependent B — )■ K*i~^i~ observables by 
the same collaboration. 

We consider the following observables: BR{B X^^y), Sk-j, Ai{B K*'y), BR{B 
Xgfi'^fi'), (v4fb)[i,6] and (-Fl)[i,6]- The experimental situation is summarised in Table IH 
together with the SM predictions. As discussed in Refs. [SI E], these observables can be 
classified as class-I (dependence only on Q ), class-II (dependence only on Cj ,Cg ,CIq) 
and class-Ill (depending on all these plus other operators, e.g., scalar operators). The 
analysis is divided into four NP scenarios: 

• Scenario A: New Physics in C7 and Cj only, real values only. 
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Figure 4: 68.3% and 95.5% CL constraints on 5Cj{^i,), 6Cj{fib) from class-I observables 
(left) and Scenario A including all observables (class-I and class-Ill, right). The con- 
straints correspond to BR{B X^'y) (purple), SK*'y (green) and Aj{B — )■ K*'y) (yellow) 
in the left panel, and BR{B — )■ X^fi^fi') (yellow), (^fb)[i,6] (orange) and (Fl)[i,6] (gray) 
in the right panel. The cross indicates the position of the benchmark point a. The origin 
of coordinates (0, 0) corresponds to the SM point. 

• Scenario B: New Physics in Cj, Cj, Cg, Ciq only, real values only. 

• Scenario B': New Physics in Cy, Cj, Cg, C[q only, real values only. 

• Scenario C: New Physics in Cj, Cj, Cg , C{q, real values only. 

The constraints from BR{B X,7), Sk-^, Aj{B K*-f) and BR{B X,/i+/i-) 
are implemented using the formulas presented in Sections 2.3 and 2.5 of Ref. [Hj- The 
constraints are set on the shift of the Wilson coefficients with respect to their SM value at 
the hadronic scale n^,. Concerning (Afb)[i,6] and (Fl)[i,6], "we provide the corresponding 
coefficients for the semi-numerical expressions of the integrated observables in Appendix 
|B| with an update of (y4FB)[i,6] and (Fl)[i,6] following the definition of the binned quantities 
described in Section [2J 

Inside the framework defined in Ref. the three observables BR{B — )• ^^7), SK*'y 
and Ai{B — >■ K*^) (class-I observables) are insensitive to New Physics contributions to 
Wilson coefficients other than the electromagnetic operators Cj . Therefore, the con- 
straints from these (class-I) observables are common to all NP scenarios and affect only 
Cj and Cj. The joint 68.3% and 95.5% CL constraints on C7 and Cj are shown in the 
left panel of Fig. HI together with the individual constraint from each observable. We 
find that the isospin asymmetry in K*^ disfavours the "ffipped-sign" solution for C7, and 
helps reducing the significance of the regions where \Cj\ ~ Cf^^Cy ~ 0. We recall that 
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our scenarios assume the coefficients of the chromomagnetic operators Og to be equal to 
their SM values. 

We then consider the impact of adding the class-Ill observables {BR{B Xgfi'^fi"), 
(Afb)[i,6] and (-Fl)[i,6]) in the different NP scenarios in turn. Within Scenario A, we obtain 
the joint 68.3% and 95.5% CL constraints in C7-C7 plane shown in the right panel of Fig. HI 
In this case, the B — )■ K*fi^fi~ forward-backward asymmetry disfavours strongly the two 
regions with large 60^ allowed at the 95.5% CL by class-I observables. This plot features 
also the benchmark point a, defined in Tableland used in the next section to study the 
power of g^-dependent B — )■ K*i~^i~ observables to discriminate among NP scenarios. 

The joint constraints within Scenario B are shown in Fig. |5l Also shown is the con- 
straint from Bg — )■ which is a direct constraint on Ciq. We see that at the 95.5% 
CL, there are four allowed regions in space of Wilson coefficients. The four benchmark 
points 61-64 are also indicated, one in each of the four allowed regions. Finally, in Figs. [6] 
and [7] we show the constraints within Scenarios B' and C, together with the corresponding 
benchmark points b' and c'. The Bg — ?■ fi^fi^ branching ratio constrains C[q directly in 
Scenario B', but only |Cio — CJg| in Scenario C. We plot the constraint based on Ref. [12], 
including the ©(AP^) correction needed to connect theory and experiment [191 1201 l2T] . 
(Using the value from Ref. [17] would result in a slightly tighter constraint. Since there is 
only an experimental bound on this branching ratio, and since the status of the theoretical 
prediction is unclear, we refrain from including this piece of information in our combined 
constraints on the Wilson coefficients.) 

As a final comment, we note that the "flipped-sign" solution for C7 is in general 
disfavoured, but not very significantly. In fact, in Scenario B, due to the loosening of the 
constraint from (Afb)[i,6] when Cg, Ciq are allowed not to vanish, this flipped-sign solution 
reappears inside the 95.5% CL. region. More precise constraints from Aj{B — t- K*'~f) and 
{Afb)ii,6] should help to settle this question. 

5 Complementarity of observables for NP studies 

The complementarity of the different angular observables in the identification of possible 
NP effects is manifest even if we deal exclusively with observables integrated over the whole 
[1, 6] GeV^ ranges. For example. New Physics contributing predominantly to Ciq will most 
likely push substantially (-P4)[i,6] below its SM valu^ (see Ref. [30]). However, this will be 
essentially indistinguishable from a NP contribution to C[q, which has a very similar effect 
on P4. This ambiguity can be resolved looking at the measured value of (-Pi) [1,6], because 
a New Physics contribution to C[q can enhance this observable considerably above its SM 
prediction, while Ciq has a negligible effect. 

Another example would correspond to a moderate enhancement of (-Pi)[i.6]- Assuming 

^To be specific, since this discussion is for illustrative purposes only, we focus on the New Physics 
scenarios of Ref. [30]. Similar conclusions can be reached by looking at Figs. [8] and |9] in this section. 
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^7 £C, SC, 

Figure 5: 68.3% and 95.5% CL contours for 5C-j{jjib)^ ^CjifJ'b), SCg^fib), CioS{fib) in Scenario 
B. The crosses correspond to the benchmark points 61 (green), 62 (blue), 63 (purple) and 
64 (orange). The blue band corresponds to the Bg — ?■ constraint. 

no significant deviation is seen in (P4)[i,6], this could signal a non-SM value of Cg or C7, Cj 
(in the island around the SM point in the C7-C7 plane). The former case will not give 
any substantial deviation in (P5)[i,6], while the latter tends to increase (P5)[i,6] above its 
SM value. Comparing (Pi)[i^6], (-P4)[i,6] and (P5)[i,6] would thus help to distinguish among 
these scenarios. 

The discriminating power of the observables is substantially increased when we con- 
sider, not only different angular observables, but also their g^-dependence. A New Physics 
contribution to C[q (in the scenario considered above) will generally increase (Pi) [1,6], but 
its effect on (P4)[i,6] and (P5)[i,6] could be washed out by simultaneous contributions to Cg 
and CiQ. However, if this enhancement of (Pi) [1,6] is also present in the region ~ 4 — 6 GeV^ 
(for instance with (Pi) [4.3,6] shifted upwards with respect to the SM), the NP effect cannot 
be misinterpreted as a non-SM value of Cj^ or Cg, because these would only enhance Pi 
in the region ~ 1 — 3 GeV^. 

These considerations apply equally well in terms of constraints when no deviation from 
the SM is observed. In Section Owe will see explicitly how different observables constrain 
the NP contributions to the Wilson coefficients. It is already useful to build an intuition 
on how different observables in different g^-bins are affected by shifts in the different 
Wilson coefficients, in order to have a better idea of the most promising observables in 
each case. Here we study briefly the effect of different New Physics scenarios on the 
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^7 ^^'J 

Figure 6: 68.3% and 95.5% CL contours for 5C-j{jjib)^ ^CjifJ'b), ^C'dif^b), ^C'^Q^fib) in Scenario 
B'. The cross is benchmark point b'. The blue band corresponds to the Bs — ?► fJ^^fJ^^ 
constraint. 

integrated observables studied in Section [31 We focus on a set of "benchmark points" 
consistent with BR{B X^^y), Sk*^, Aj{B K*^), BR{B Xsfi'^iJ,'), (v4fb)[i,6] and 
{Fl)[i,6], according to the analysis of Section HI These points are specified in Table El and 
also indicated in Figs. HHIl 

Each NP benchmark point gives a prediction for each observable as a function of g^, 
shown in Figs. El and [HI together with the SM predictions. In these figures, the yellow and 



Point 




SC'-^inb) 


SCg{flb) 


5Cio(/ife) 


SCgiHb) 




a 


-0.04 


-0.1 














bl 


-0.03 


-0.09 


-1.5 


-1 








62 


0.3 


-0.4 


1 


6 








63 


0.45 


0.45 


-9 


2 








M 


0.9 


0.1 


-9 


8 








b' 


-0.05 


-0.15 








3 


1.5 


c 


0.4 


-0.45 


-7 


5 


3 


5 



Table 5: NP benchmark points used in the analysis of Section O 
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Figure 7: 68.3% and 95.5% CL contours for SC7{iJ,b), SC'-^{iib), SCg {fj^b), 5C}o (^^b) in Scenario 
C. The cross is benchmark point c. The blue band corresponds to the Bs — t- fi^fi^ 
constraint. 
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Figure 8: Binned SM predictions for the observables (Pj) together with NP curves corre- 
sponding to the benchmark points a (blue), bl (green), 62 (red), 63 (purple dashed), 64 
(orange dotted), 6' (black dot-dashed) and c (black dotted), given in Table [51 
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Figure 9: Binned SM predictions for the observables (Afb) and {Fl) together with NP 
curves corresponding to the benchmark points a (blue), 61 (green), 62 (red), 63 (purple 
dashed) , 64 (orange dotted) , 6' (black dot-dashed) and c (black dotted) , given in Table 



red rectangles correspond to the SM value, like in Figs. [H |2]and[3l The curves correspond 
to the NP scenarios represented by each benchmark point collected in Table [51 We can 
draw a few conclusions: 

• The points 62, 63 and 6' enhance the observables (Pi) and (P4) substantially, and 
the impact on (Pi) ((P4)) is more important for > 3 GeV^ (g^ < 3GeV^). 

• (Pi) and (P4) do not allow one to discriminate among 62, 63 and 6', but (P2) and 
(P5) exhibit distinctive behaviours for 62 and 63. For instance, a suppression of (Pi) 
below ~ —0.5 together with an enhancement of (P2)[4,6] above ~ would favour 
62 and 63 with respect to 6'. This conclusion could be verified by a suppression of 
(P2)[i,3] below its SM value and an enhancement of (P5)[4,6]- 

• A similar situation occurs with the points 64 and c. The observable (P2)[5,6] could 
favour these scenarios, but cannot distinguish among them. However, (Pi) [1,2] and 
(Pi) [5 6] can discriminate these scenarios, as well as (Pg)[i_6] if the experimental values 
are accurate enough. 



6 The benefit of using clean observables 

In this section we discuss the advantages of using clean observables in analyses of P — )► 
K*i~^i~ as opposed to other observables such as 5*3 or Aim. For definiteness we focus on 



23 




(GeV-) if (ceV-) 

Figure 10: Predictions in the SM and in the case of NP at the benchmark point 62 for Pi 
(left) and 6*3 (right). The yellow boxes are the SM predictions integrated in five 1 GeV^ 
bins. The blue curve corresponds to the central values for the NP scenario. The green 
band is the total uncertainty considering the form factors of Refs. [261 [23], while the gray 
band is the total uncertainty obtained using the form factors of Ref. |25]. In the case of 
Pi the gray band is barely visible. 



Pi and 53, but it should be kept in mind that the conclusions are more general. 

As discussed extensively in Refs. [281 ISH EH [30], clean observables are constructed 
in such a way that at LO and at large recoil, an exact cancellation of the form-factor 
dependence occurs. This indicates that clean observables should be stable under variation 
of hadronic uncertainties, as opposed to other observables, such as F^, ApB, S3, etc. This 
is relevant because of the spread of published errors in the determination of form factors 
from light-cone sum rules (see Refs. [251 [26] and the introduction). 

If the form factors of Ref. [25] are used in the evaluation of Fl, for example, the error 
bars get enlarged by a factor of three (this occurs also for other observables with a strong 
sensitivity on form factors, like Afb)- On the contrary, this enlargement does not happen 
in the case of Pi, which is practically insensitive to these uncertainties. In the case of 5*3, 
an accidental circumstance makes its SM uncertainty smaller than what one would infer 
from the fact that 5*3 ~ PiFt (that is, a similar percentual enhancement of the errors as 
Fl). The fact that Pi ~ in the SM, makes 5*3 almost insensitive to the error in Fl only 
near the SM point. This makes the situation with a bit more subtle. The important 
point here is that in the presence of New Physics, an enhancement of Pi produces an 
enlargement of the error bars in the theoretical prediction for 5*3 automatically, which 
makes 5*3 almost unable of discriminating between NP models where Pi does not vanish. 

In Fig. [To] we show the SM predictions and the predictions for benchmark point 62 
(see Table [5]) for Pi and S3 calculated with both choices of form factors (Refs. [25] and 
[26]). We find that: 
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• The SM prediction for Pi is insensitive to the choice of form factors. 

• The SM prediction for S3 shows a mild dependence on the choice of form factors, 
and hadronic uncertainties are enlarged up to a 50% when using the form factors of 
Ref. [25] compared to those in Ref. [26] . 

• The NP prediction for Pi is insensitive to the choice of form factors. 

• The NP prediction for 5*3 is very sensitive to the choice of form factors. Indeed, the 
hadronic uncertainties increase from a factor 2 to a factor 3 when using the form 
factors of Ref. [25] . 

The conclusion is that, in New Physics analyses of 5 — > K*i^i^, when using the ob- 
servable 5*3, one is obliged to take into account hadronic uncertainties at each point in 
the NP parameter space, and these uncertainties can be substantial. On the other hand, 
the observable Pi is not affected by this disease, and it is therefore a much more robust 
observable. Moreover, from Fig. [10] we can see that Pi has the potentiality of completely 
excluding a NP scenario given by benchmark point 62, while S3 can barely distinguish 
this point from the SM case at more than 2a. We expect the same results to hold in other 
regions of the NP parameter space. 

We would like to emphasise that similar conclusions are expected for other form factor 
dependent observables. For example, the observable Aim ~ 5*9 ~ PjPt will also be 
protected from large hadronic uncertainties near the SM point, since P3 ~ in the SM. 
However, complex NP can enhance considerably P3 while being consistent with other data 
(see Ref. [30]). In these NP scenarios, A^-^ is expected to suffer from a similar problem as 
the one described for S3 (a problem that does not affect P3). Furthermore, such problems 
may also happen also in the case of CP conjugated observables such as Ag and ^3. 

In view of this situation, one should be particularly careful with a rather usual hy- 
pothesis consisting in assigning the same hadronic uncertainty on observables in any NP 
scenario as in the SM case. This might be a good approximation in the case of clean 
observables (see for instance the left-hand side of Fig. [TOl where the uncertainty on Pi 
is indeed similar from the SM case to the NP scenario considered), but it can be very 
misleading for other quantities, sensitive to form factors (as illustrated on the right-hand 
side of Fig. [TOl where the uncertainty on S3 is significantly enlarged from the SM value 
to the NP scenario considered). 

7 Extracting clean observables from existing experi- 
mental measurements 

As pointed out in Section [6l there is a clear advantage in using the clean observables 
-Pi, 2,3 rather than 5*3, Aim and Ap^ (or F^) to put constraints on the Wilson coefficients. 
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n ««+0.06+0.04 
"-'■"J"J-0.06-0.03 




n no+o.io+o.os 

"^•"^^-0.06-0.01 


n 09+0-07+0.01 
"-'■"-'^-0.07-0.01 


n 07+0.07+0.02 
"-'■"-" -0.07-0.01 




n nc^+o-18+0.05 

'-'•'-'^-0.12-0.01 


n -I O+0.13+0.01 

"J- -"-0-0.13-0.01 


n in+o.i5+o.o2 

"-'■-'-"-'-0.16-0.01 



Table 6: LHCb experimental results for binned observables (from Ref. |4j). 



However, this can be achieved if the experimental fits are performed consistently consid- 
ering these observables (see Ref. [30]). In its latest experimental analysis of g^-dependent 
observables in i? — )■ K*£^£~ , the LHCb collaboration [3] has preferred to fit directly the 
coefficients of the angular distributions, providing only observables proportional to the Jj 
coefficients (in particular (5*3) and (Ai^)). We have collected for reference the relevant 
experimental results for these observables in Table [61 

Of course, we can compute the clean observables from the measurements provided 
using the formulas of Sections 3 and 4 of Ref. [30], or equivalently from Eqs. (|5]). If this is 
done without knowing the correlation matrix (which is not provided by the experimental 
collaborations yet), one obtains errors much larger that the real uncertainties. It is still 
worth exploring the current situation based on these observables while waiting for cor- 
related values. In addition, higher-statistics analyses from LHCb are expected to reduce 
the experimental errors on these observables considerably very soon. 

Attending to these considerations, experimental values for -Pi, 2,3 can be derived from 
the measurements of 6*3, A^^, Ap^ and Fl in Table [6] by means of Eq. ([5]): 

/P\ _ 2(5*3)^11 ,p> _ 2 (v4FB)bin ,p> _ (Am) bin 

^ " 1 - (FL)bin ' ^ " 3 (1 - (F,)bin) ' ^ '^'^'^ " (1 - (FL)bin) ' ^ ' 

In Table[7|we present the resulting experimental values for (Pi), (-P2), (-P3) in the different 
bins, together with their SM predictions. The experimental errors are calculated in the 
following way. We first add in quadrature both errors in Table [6l and symmetrise upper 
and lower uncertainties. Assuming these errors are Gaussian, the errors for -Pi,2,3 are 
obtained by the usual error propagation formula from Eqs. ([2j 
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Observable 


Experiment 


SM prediction 


(-Pi) [2,4.3] 


-0.19 ±0.58 


-0.051 ±0.050 


(-Pi) [4.3,8.68] 


0.42 ±0.31 


-0.117 ±0.059 


(A) [1,6] 


0.29 ± 0.47 


-0.055 ± 0.051 


(-P2) [2,4.3] 


0.51 ±0.27 


0.232 ±0.069 


(-P2) [4.3,8.68] 


-0.25 ±0.08 


-0.405 ± 0.064 


(^2) [1,6] 


0.-35 ± 0.14 


n no A 1 n c\nn 

0.084 ± 0.066 


(-P3) [2,4.3] 


0.08 ±0.35 


-0.004 ± 0.024 


(-P3) [4.3,8.68] 


-0.05 ±0.16 


-0.001 ±0.027 


(^3) [1,6] 


-0.21 ±0.21 


-0.003 ± 0.024 



Table 7: Experimental values for the clean observables Pi, Pi and P3 within different g^- 
bins, extracted from the measurements of 5*3 , Aim, ^fb and Fl, and their SM predictions. 

8 Present and future constraints from Q'^-dependent 
B K*i^i~ observables 

The recent LHCb measurements for g^-dependent B — )■ K*i^i~ observables ^ are divided 
into four bins -if we restrict ourselves to the low-g^ region. These bins are [0.05, 2], [2, 4.3], 
[4.3,8.68] and [1,6] GeV^. These results yield a first glimpse of the future, where precise 
measurements of the full angular distribution within fine g^-bins will be available. The 
purpose of this section is to study the impact of the g^-dependent observables provided in 
Ref. [1] on the constraints on the Wilson coefficients, and to analyze what are the future 
expectations concerning the constraints from these observables. 

A brief discussion is in order concerning the results in Ref. [4] . From the theory point 
of view the first bin [0.05, 2] is very difficult to control, since the decay rate contains 
contributions from light resonances below g^ ~ IGeV^. The third bin [4.3,8.68] is also 
more difficult to handle theoretically, as it gets near to the cc resonance region. For 
this first attempt, we choose to drop the first bin and to consider constraints from the 
two others [2,4.3], [4.3,8.68]. We do not include [1,6] as the results are likely to be 
strongly correlated with the two smaller bins considered for this study (we remind that 
the averaged experimental results in the bin [1, 6] have already been considered in Section 
Hlin the case of {Ap^) and (Pl))- 

8.1 Constraints from binned (^fb) and (F^) 

We first consider the observables (Afb) [2,4.30] , (^fb) [4.30,8.68], (P'l) [2,4.30] and (Pl) [4.30,8.68]- 
The experimental numbers for these observables are given in Table El The constraints are 
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6Cj ^Ct 



Figure 11: Left: Individual 68.3% and 95.5% CL constraints in the 5C7{fib)-SCj{fib) plane 
from (Afb) [2,4.30] and (Afb) [4.30,8.68], together with the combined constraint. Right: Same 
analysis for (Fl) [2,4.30] and (i^i) [4.30,8.68]- The origin corresponds to the SM value. 

implemented using Eq. in App. IB] together with the coefficients in Tables |8] and [91 

The individual 68.3% and 95.5%C.L. contours for these observables in the C7-C7 plane 
are shown in Fig. [11] We see that: 

• The constraints from (Afb) are consistent with the SM at 95.5%C.L. Some tension is 
caused by (^fb) [2,4.30], while (^fb) [4.30,8.68] is compatible with the SM at 68.3%C.L. 
This tension disappears if the form factors of Ref. [25] are used in the SM predictions. 

• The constraints from (Fl) show a discrepancy with the SM, which is outside the 
95.5%C.L. region. Both observables (Fl) [2,4.30] and (Fl) [4.30,8.68] have a tendency to 
avoid the SM point, although this feature is stronger in the case of (Fl) [4.30,8.68]- 
Again, this tension goes away if the form factors of Ref. [25] are used for the SM 
predictions. 

• Taking seriously these tensions would require a NP that affects simultaneously Aps 
around ~ 3 GeV^ and Fl around ~ 6 GeV^. 

8.2 Constraints from binned (Pi), (P2) and (P3) 

We now consider the constraints from the observables (Fj) [2,4.30] , (-Pj) [4.30,8.68], {Pi) [1,6] 
with i = 1,2,3. The experimental values and SM predictions for these observables are 
collected in Table [71 The constraints are implemented using Eq. f[M|) in Appendix [B] 
together with the coefficients in Tables [TOl [HI and [121 
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Figure 12: Left: Individual 68.3% CL constraints in the SCi^jibj-SCj^iJih) plane from the 
integrated clean observables (A) [2,4.30] and (Pi) [4.30,8.68], together with the combined result. 
The red region and contour correspond to the combined 68.3% and 95.5% CL regions. 
Right: Same analysis for (P2)bin- The origin corresponds to the SM point. Currently, 
(-P3)bin does not provide any constraint on the SC7{fib)-SC'-^{fib)- 

The individual constraints from these observables in the C7-C7 plane are shown in 
Fig. [121 We see that: 

• The constraints from (Pi) are not very stringent yet. However there is a mild dis- 
crepancy of (Pi) [4.30,8.68] with respect to the SM, as well as the combined constraint 
from the two bins. This result is not affected by form factor uncertainties. 

• The constraints from (P2) are already quite interesting. The two bins point towards 
a negative value of 6C7, and the SM point is just in the border of the 95.5% CL 
region. Again, this result is not affected by form factor uncertainties. While the 
theoretical prediction for (P2) [4.30,8.68] can suffer from the proximity of the bin to the 
cc resonance, we point out that the same tendency to negative 6Cy is hinted at by 
the observable (P2)[i,6], indicating that this is not a feature introduced by the data 
above 6 GeV^. 

• The constraints from (P3) are inconclusive for the time being. This could be guessed 
already from the NP analysis of Section O (see Fig. |H]). It is well known that the 
CP-averaged version of (P3)bin (the one we are considering here) is not very sensitive 
to NP, and that the corresponding CP-asymmetry might be more interesting when 
constraining NP (see for example Ref. [12] )■ In this case we would suggest to focus, 
instead of Ag, on the corresponding clean CP-asymmetry, since Ag can be affected 
by the problems discussed in Section [6] concerning its sensitivity to form factors. 
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-OS 0.0 0.5 1.0 1.5 

Figure 13: Individual constraints in the C7 — Cj plane from hypothetical measurements 
of the observables (Pi)[2,4.3] (yellow), (P2)[2,4.3] (green), (Pi)[2,4.3] (purple) and (P5)[2,4.3] 
(blue), corresponding to central values equal to the SM predictions and an experimental 
uncertainty aexp = 0.10. The combined 68.3% (dark red) and 95.5% (light red) C.L. 
regions are also shown. 

We stress again that these constraints should be considered as conservative, since they 
are based on the experimental numbers extracted in Section [7| in absence of experimental 
correlations. Therefore, the uncertainties of Pi, 2.3 that we quote in Table [7] are most 
probably overestimated, and could be reduced significantly once experimental correlations 
are available. 

8.3 Future prospects 

The experimental numbers in Table [6] for the various observables, as measured by the 
LHCb collaboration, contain uncertainties at the level of ~ 0.10, and up to ~ 0.20 for 
the observable {S3). The numbers for (Pi, 2,3) extracted in the previous section from the 
measurements in Table [6] contain larger uncertainties, up to ~ 0.5 (see Table [7j). As dis- 
cussed above, these errors are probably overestimated since they do not take into account 
the relevant correlations among the observables. It is reasonable to expect that a direct 
extraction of (Pi, 2,3) from the data would give, with the present statistics, experimental 
uncertainties for these observables in the ballpark ~ 0.10 — 0.20, as is the case for the 
observables in Table El 

We have also seen in Section 18.21 that even with the enlarged uncertainties of Table Ui 
the constraints on C7, Cj from (P2) are nontrivial. These uncertainties will improve very 
significantly in the near future in forthcoming analyses of larger data sets by the LHCb 
collaboration. 
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In order to illustrate the very large impact that moderately precise measurements of 
the clean observables will have on New Physics, we consider the constraints on C7, Cj 
by hypothetical measurements of (-Pi) [2,4.3], (-P2) [2,4.3], (-P4) [2,4.3] and (P5) [2,4.3] with central 
values at their SM points, and experimental uncertainties of aexp = 0.10. It is important 
to emphasise that these errors are not only realistic but also achievable with the current 
statistics. The result of this exercise is presented in Fig. [131 where the individual con- 
straints in the C7 — Cj plane (corresponding to Scenario A) from these four observables 
are shown together with the combined 68.3% and 95.5% CL. regions. Clearly these ob- 
servables will play a very important role in the future, with the potentiality of providing 
the first unambiguous signal of New Physics in the flavour sector. 

9 Comparison with other works 

Several theoretical studies [HI O El IH [TOl [TTl [121 150] have exploited the recent measurements 
of 6 — )• 57 and b — )■ si~^i~ . We mention here briefly some differences of recent references 
with respect to the present work where only low-g^ bins were considered. With that 
respect, one should notice that these references use at most the B — )■ K*i'^i~ observables 
integrated over [1,6] GeV^, but not the values on finer bins included here. 

In Ref. [50], a detailed discussion was provided concerning additional observables re- 
lated to the photon polarisation in 6 — > 57 transitions, allowing one to constrain C7 
and C7 more precisely. The three processes of interest were — )■ i^"*°(— )■ Ks'n'^)'y, 
B -ft'i(-)- Ktx'k)'-^ and 5° K*^{-^ K^T[~)l'^i~ at low g^, showing the potential 
of an analysis combining all four processes even in the case of complex Wilson coeffi- 
cients (assuming that there are only small contributions from the other operators for 
B^ — )■ K*°(— )■ K^Ti^)l^ l'^ in this energy region). In particular, the current constraints 
on CyC) from B{B — )► Xs^i) and S'x*7 were studied, either with real or complex values 
for these Wilson coefficients (and SM values for the other operators). In the real case, a 
four-fold degeneracy is observed, corresponding essentially to the regions obtained in our 
Scenario A without B — )■ i^*7 isospin asymmetry (see fig. [1]). These ambiguities (which 
are even more numerous in the complex- valued cases) can in principle be lifted once more 
observables are measured from all three processes, with different strengths and weaknesses 
depending on the NP scenario considered. 

In Ref. [12] (which updates Ref. [9]), branching ratios for the inclusive modes B — ^ 
Xs7 and B — t- Xsd-d, (low and high-g^) and the inclusive CP asymmetry for b ^ s'j are 
considered, as well as the branching ratio for B — )■ Kii and Bs — )■ fi'^fi" and observables 
for B — )■ K*ii (branching ratio, longitudinal polarization fraction, forward-backward 
asymmetry, Ag and 5*3 both in low- and high-g^ regions). The constraints are put on 
C7 ('), CqC), Cio('), CsC), Cp(') with real or complex values, first in a frequentist framework, 
then in a Bayesian approach, with a focus on some specific NP scenarios. Once again, a 
good agreement with the SM is obtained. In addition, the possibility of a sign-flip in C7 
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alone is disfavoured due to the branching ratio of S — Xs7 as well a.s ApB- In the scenario 
without right-handed currents and with real values of the coefficients of the SM operators, 
the flipped sign solution Cy g^io — > —^7,9,10 remains allowed (since the B — )■ K*'y isospin 
asymmetry was not included in the analysis). In this scenario, the constraints obtained on 
the Wilson coefficients are similar to those that we obtain for scenario B, up to additional 
solutions with C7 7^ 0. The very large parameter space and the choice of different scenarios 
prevents us from comparing our results in more detail, but we stress that the analysis in 
Ref. [12] includes form-factor sensitive observables like 5*3. 

In Ref. [11], the authors did not consider the inclusive modes B — )■ Xg'-f and B — t- Xgii. 
On the other hand, they considered the branching ratio of the exclusive modes B — )■ K*^, 
B — )■ Kii, Bs — fi^fi~ as well as the available observables for B — )■ K*'y (branching 
ratio, longitudinal fraction, forward-backward asymmetry, transverse asymmetry Aj, and 
5*3 both in low- and high-g^ regions). The constraints were analysed in the SM basis, 
constraining only real values of Cj,Cg,CiQ, in a Bayesian framework. The inputs for the 
hadronic form factors are taken from the light-cone sum rule analysis of Ref. [26]. They 
provide also predictions for the other transverse asymmetries A^j) , as well as for the low- 
recoil observables H^^ introduced in Ref. p3]. Obviously, as the list of inputs, scenarios 
and statistical frameworks are rather different, we can only perform a limited comparison 
with our benchmark points in different NP scenarios. In both analyses, the SM point is 
favoured. A second solution, with flipped signs Cy^g^io — )■ — Cy^g^io is allowed in Ref. [11] 
as in the previous references. The authors observe an update from prior to posterior 
p.d.f.'s of the form factors, which can be interpreted as the fact that the data themselves 
constrain partly the hadronic uncertainties, due to a slight tension between B — )■ Kii 
and B — t- K*ii. In view of our discussion concerning theoretically clean and form-factor 
sensitive observables, it would be interesting to perform a similar analysis separating the 
two sets of observables. 

10 Conclusions 

The decay mode B — )■ K*i^i~, with its large set of angular observables, is becoming a 
more and more important process in constraining New Physics models. These constraints 
are complementary to those from inclusive and exclusive radiative decays. With increas- 
ingly precise experimental data on these modes and the prospects for the future from the 
LHC, theoretical uncertainties must be kept under control. In this context, the use of 
theoretically clean observables is not only convenient but also extremely recommendable, 
or even mandatory. The sensitivity of different observables to hadronic uncertainties has 
been addressed in Section |6l and the conclusion is that clean observables such as Pi are far 
more robust than other observables like Fl or S3, translating into a better performance 
in discriminating among different models. 

A full description of the angular distribution of the B — )■ K*i^i~ decay in terms of 
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a maximum set of clean observables was presented in Ref. [30], where the observables 
-Pi, 2,3,4,5,6 where introduced. In this paper we have presented a simple and compact ex- 
pression for the coefficients of the distribution in terms of these clean observables, and 
given SM predictions for these observables integrated over a series of g^-bins, that can 
be directly compared with experimental data. These predictions are collected in Table 
[1] and in Figs. [1] and [2j As an important point, we have seen that the first three clean 
observables -Pi, 2,3 are already measured, and can be extracted form the latest measure- 
ments by the LHCb collaboration [4J. The experimental numbers for these observables 
together with their SM predictions are given in Table [71 The uncertainties attached to 
these numbers are certainly overestimated, as we did not have the required correlations 
among experimental measurements (and thus treated them as uncorrelated) . 

We have also studied the model-independent constraints on the Wilson coefficients 
Q , Cg , from radiative decays B — )■ Xg'y and B — )■ K*'j, and semileptonic decays 
B — !■ Xsfj,^fj,~ and B K*fi~^fi~. Excluding all B — )• K*fi~^fi~ observables except for 
the integrated observables (Afb)[i,6] and (-Pl)[i_6] leads to the set of constraints shown 
in Figs. miZl Identifying a set of benchmark New Physics points compatible with those 
bounds, one can see that very large New Physics contributions to other observables in 
B — 7- K*fi~^fi~ are allowed, specially in Pi, P2, -P4 and Pg (see Fig. [8]). We have investigated 
the present constraints imposed by Pi, P2 and P3 on C7, Cy, with already quite interesting 
constraints from P2 suggesting a lower value of Cj than the SM value. We then showed 
the powerful prospects that the set of clean observables Pi,2,3 and Pi'^s'^e' will have in the 
short term to discriminate possible New Physics contributions, illustrated in Fig. [131 

Considering the advantages provided by the use of the Pj observables at large recoil, 
we hope that the present study will be a strong incentive for experimentalists to rephrase 
their study of the low-g^ region in terms of these observables. 
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A Jg in terms of Pi observables and the Q observable 

In Section[2lwe have provided the explicit expressions for the coefficients of the distribution 
in terms of the observables of the basis, neglecting mass terms and scalar contributions|f| 
In this case, one finds very simple and compact expressions for all these coefficients (see 
Eqs. flT^ ) with the exception of Js- The reason is that Jg, in the absence of scalar con- 

■^The general case including lepton masses and scalars is discussed in full detail in Ref. [3D]. 
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tributions, is not an independent quantity (exactly like Ju = 3J2s and Jic = — J2c in the 
massless case) and deserves a separate discussion. The counting of degrees of freedom 
and continuous symmetries in this case shows that there are only 8 degrees of freedom 
parametrised by the observables Pi=i„s, and dV/dq^. This means that Jg can be 
expressed in terms of these observables: 



^^^|(P2P6 - P3P4) + V - + ^5(^2^4 + P3Pe)^Jl-P? 



Jr — —\ 



+ \a-t./"^(^-p"^-p"-p"y] (28) 

This expression is found by solving Jg in terms of the other coefficients using the relation 
in Eq. (3.15) of Ref. [31] together with Eqs. ( !T4|) . One can also replace -P4,5,6 by -Pi 55 
using Eqs. ffT5|) -f lT71l . 

We notice that in Eq. ( l28ll . a discrete quantity rj is left as a free parameter that can 
take only two values ±1 in the massless cas^. This parameter is indeed an observable, 
and its SM prediction is rj^'^ = +1 as can be seen by substituting the SM values for the 
observables Pi. This has the interesting consequence that a measurement of7] = —1 would 
be an unambiguous indication of New Physics (originating, for instance, from new weak 
phases or sign flips in Wilson coefficients). Deviations from |?7| = 1 could be expected, 
for instance, from scalar contributions entering P5 and Pq, since such contributions would 
break the symmetry relation among the coefficients of the angular distribution. 

However, when one tries to write a similar relationship for "binned" observables, it is 
clear that a naive substitution Pi — )■ (Pi)bm is not possible due to the highly non-linear 
form of Eq. ( |28l) . One simple solution is to introduce an extra clean observable Q: 

Js = -\Q'^F^l^2 ' (29) 
with Q' = (5vT+~Pi" and whose definition, also valid in the massive case, is 

^ Im(n^nx) ^ V2J8 

The vectors nj (with i = 0, _L) are defined in Ref. [3^. This clean observable Q (or Q') is 
related to the form factor-sensitive observable Ss by means of 

Q' = . (31) 



"'if one solves Eq. (3.15) of Ref. [28] for any coefficient of the distribution whose global sign could 
depend on the model (SM or NP), like J4, J5, Jg^, J7, Jg or Jg, it should yield a free parameter rj for 
each J (which are not independent of each other). 
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Figure 14: Standard Model prediction for the observable Q', unbinned (continuous curve) 
and binned according to 5 bins of width 1 GeV^. 

For completeness we also provide the expression of this observable in terms of our basis 
of observables: 



Q 



{P^Pq - PsPa) + V (^2^6 - PsPiY + ^5(^2^4 + P3Pe)\n - Pi 



(32) 



In principle, this implies that Q (or Q') can be measured using either Eq. ( !30!) or Eq. (!32l) 
(with 7] = +1 in the SM but free in general). However, Eq. fl32|) would be of practical 
experimental use only in the limit of the size of the binning going to zero. Of course, this 
expression Eq. fl32|) can be very easily turned into an equation that relates 5*8 with the 
other Si by using Eq. (|5]) and Eq. (fT8|) . This shows the redundancy of in the massless 
case, again up to a single discrete parameter rj. 

At a more practical level, the integrated form of this observable from Eq. ( 130|) is: 



{Q') 



bin 



bin 



\//bin^^^C4(g^) /bi^c/g2(co(g2) _ ^^(^2)) ^(C4)bin ( (Co)bin " (C4)bin) 



(33) 



and its SM prediction is given in Fig. [TH Notice that, as mentioned above, the binning 
procedure breaks the relation Eq. (3.15) of Ref. [31] among the coefficients of the distri- 
bution, which is only recovered in the limit of the size of the binning going to zero where 
the observable rj would become accessible. 
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In conclusion the set of Eqs. (fT^ and Eq. ( p8|) provides a complete parametrisation of 
the distribution in the massless case with no scalars in terms of only eight parameters 
(eight observables -Pi, 2,3,4,5,6 and Fi^dV /dq^) in agreement with the eight independent 
degrees of freedom, but there is one extra redundant Q (or Q') observable (up to a 
discrete parameter rf) which however can be fixed only once (J8)bin is measured, due to 
the binning procedure adopted by experimental analyses. 



B New Physics expressions for binned and unbinned 
observables 

In this appendix we present the numerical expressions for the integrated observables 
(Afb), {Fl), (-Pi, 2,3) and (-P^se) as a function of the NP Wilson coefficients, for different 
choices of the g^-binning. 

The formula for any observable {Ok) has the general form: 

\Ok) = ^ — Jr~^ ^ ^ > 

Z^i,jr=0,7,7',9,9', 10,10' -^(i j) "'"'i 

where 5Cq = 1. The coefficients 6, Nij and Dij corresponding to each observable are 
collected in Tables |81IT2| where only nonzero coefficients are displayed. The coefficient 6Ci 
denotes the NP contribution to the Wilson coefficient Ci at the hadronic scale fib- The 
parameters 6k correspond to the theoretical error assigned to each observable, and are 
also collected in Tables [8lfT2l We assume here that the uncertainty obtained within the 
SM can be considered as a good estimate of the uncertainty for arbitrary (NP) values of 
the Wilson coefficients. As discussed in detail in Section [6l this assumption is appropriate 
for theoretically clean observables only. 



C Statistical approach 

We determine and combine our constraints in a frequentist framework, treating theoretical 
and experimental uncertainties on the same footing (i.e., taking them as normally dis- 
tributed random variables). The model independent bounds are obtained in the following 
way. A chi-square function is constructed, according to: 

where p are the theoretical parameters constrained by the analysis (Wilson coefficients, 
hadronic quantities, CKM matrix elements. . . ), Of^, 6k are the central value for theoretical 
prediction and error for the observable Ok, and O^^^ ± ak is the experimental average. In 
case of several uncertainties (statistical and systematic ones, for instance), we combine 
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-103.17 


-107.05 


^(7', 9') 


1292.76 


592.97 


1175. 


246.66 


254.16 


259.91 


264.4 


267.64 


^(9,9') 


-220.54 


-100.27 


-214.08 


-40.28 


-42.38 


-44.26 


-46.01 


-47.61 


-'^(10,10') 


-220.54 


-100.27 


-214.08 


-40.28 


-42.38 


-44.26 


-46.01 


-47.61 


-D((),0) 


4339.85 


1872.7 


5075.88 


854.2 


790.49 


823.22 


890.38 


981.56 


D((),T) 


45.73 


152.6 


3422.74 


-1041.84 


-170.54 


211.26 


439.85 


607. 


-C{0,7') 


-1738.05 


-781.35 


-1994.7 


-294.31 


-324.1 


-348.2 


-372.01 


-399.42 


-0(0,9) 


774.47 


333.24 


1118.83 


106.28 


129.25 


153.17 


178.6 


207.17 


-D{0,9') 


-814.41 


-366.11 


-903.42 


-145.14 


-153.63 


-162.12 


-171.24 


-182.28 


-C{0,10) 


-1200.92 


-535.82 


-1323.08 


-194.54 


-218.18 


-241. 


-263.07 


-284.13 


r-i 


you.zz 






1 ( o.Oo 


ioz.oy 


1 on '71 




zUo. io 


-0(7,7) 


15231.1 


6579.08 


6569.48 


5604.07 


3420.87 


2513.2 


2008.91 


1684.05 


-C{7',7') 


15289.4 


6605.3 


6620.35 


5616.8 


3432.41 


2524.48 


2020.21 


1695.49 


-0(9,9) 


278.73 


124.36 


307.08 


45.15 


50.64 


55.93 


61.06 


65.95 


-0{9',9') 


278.73 


124.36 


307.08 


45.15 


50.64 


55.93 


61.06 


65.95 


-0(10,10) 


278.73 


124.36 


307.08 


45.15 


50.64 


55.93 


61.06 


65.95 


-0(10', 10') 


278.73 


124.36 


307.08 


45.15 


50.64 


55.93 


61.06 


65.95 


-0(7,7') 


-1195.66 


-540.43 


-1115.31 


-236.61 


-232.41 


-236.1 


-242. 


-248.55 


-0(7,9) 


1282.19 


588.17 


1164.65 


244.74 


252.13 


257.79 


262.19 


265.34 


-0(7,9') 


-465.63 


-211.33 


-457.65 


-84.32 


-89.08 


-93.42 


-97.5 


-101.32 


-0(7', 9) 


-493.35 


-224.04 


-482.85 


-89.61 


-94.53 


-98.99 


-103.17 


-107.05 


-0(7', 9') 


1292.76 


592.97 


1175. 


246.66 


254.16 


259.91 


264.4 


267.64 


-0(9,9') 


-220.54 


-100.27 


-214.08 


-40.28 


-42.38 


-44.26 


-46.01 


-47.61 


-0(10,10') 


-220.54 


-100.27 


-214.08 


-40.28 


-42.38 


-44.26 


-46.01 


-47.61 



Table 8: Coefficients for the New Physics formula of (^fb)- 
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[1,6] 


[2,4.3] 


[4.3,8.68] 


[1,2] 


[2,3] 


[3,4] 


[4,5] 


[5,6] 


s 


0.04 


0.035 


0.04 


0.059 


0.038 


0.033 


0.034 


0.037 


^(0 0) 


1754.68 


793.06 


1880.43 


318.5 


334.72 


350.22 


366.2 


385.04 




805.76 


363.46 


872.59 


145.27 


153.08 


160.66 


168.62 


178.13 


^(0 7') 


-823.42 


-371.42 


-891.71 


-148.46 


-156.43 


-164.18 


-172.31 


-182.04 




403.06 


181.92 


434.57 


72.88 


76.69 


80.38 


84.24 


88.87 




-403.06 


-181.92 


-434.57 


-72.88 


-76.69 


-80.38 


-84.24 


-88.87 




-464.07 


-211.81 


-438.02 


-86.04 


-89.95 


-93.28 


-96.18 


-98.62 


-^fO lOM 


464.07 


211.81 


438.02 


86.04 


89.95 


93.28 


96.18 


98.62 


^(T 7) 


430.42 


196.23 


409.76 


79.34 


83.19 


86.5 


89.44 


91.96 


N(^, 7/-, 


449.49 


204.92 


427.91 


82.85 


86.87 


90.34 


93.4 


96.03 


^(9 91 


107.71 


49.16 


101.66 


19.97 


20.88 


21.65 


22.32 


22.89 


-^(9', 9') 


107.71 


49.16 


101.66 


19.97 


20.88 


21.65 


22.32 


22.89 


-^(10,10) 


107.71 


49.16 


101.66 


19.97 


20.88 


21.65 


22.32 


22.89 


-^(10', 10') 


107.71 


49.16 


101.66 


19.97 


20.88 


21.65 


22.32 


22.89 


JV(7,7') 


-597.83 


-270.22 


-557.66 


-118.3 


-116.2 


-118.05 


-121. 


-124.28 


■^(7,9) 


641.09 


294.08 


582.32 


122.37 


126.07 


128.89 


131.1 


132.67 


JV{7,9') 


-232.82 


-105.66 


-228.83 


-42.16 


-44.54 


-46.71 


-48.75 


-50.66 


JV{7',9) 


-246.68 


-112.02 


-241.42 


-44.81 


-47.26 


-49.49 


-51.58 


-53.53 


I, ' ^ 


646.38 


296.49 


587.5 


123.33 


127.08 


129.95 


132.2 


133.82 


-^{9,9') 


-110.27 


-50.13 


-107.04 


-20.14 


-21.19 


-22.13 


-23.01 


-23.81 


-^(10,10') 


-110.27 


-50.13 


-107.04 


-20.14 


-21.19 


-22.13 


-23.01 


-23.81 


^(0 01 


2169.93 


936.35 


2537.94 


427.1 


395.25 


411.61 


445.19 


490.78 


D(Q 71 


22.87 


76.3 


1711.37 


-520.92 


-85.27 


105.63 


219.92 


303.5 


^fO 7M 


-869.02 


-390.67 


-997.35 


-147.15 


-162.05 


-174.1 


-186.01 


-199.71 


^fn 91 


387.24 


166.62 


559.41 


53.14 


64.62 


76.59 


89.3 


103.59 


.^fO 9M 


-407.2 


-183.06 


-451.71 


-72.57 


-76.82 


-81.06 


-85.62 


-91.14 




-600.46 


-267.91 


-661.54 


-97.27 


-109.09 


-120.5 


-131.54 


-142.07 


^fn io'l 


475.11 


216. 


461.19 


86.76 


91.29 


95.36 


99.12 


102.57 


D(7 71 


7615.55 


3289.54 


3284.74 


2802.04 


1710.43 


1256.6 


1004.45 


842.02 


^(7' 7M 


7644.7 


3302.65 


3310.18 


2808.4 


1716.21 


1262.24 


1010.11 


847.74 


-C{9,9) 


139.37 


62.18 


153.54 


22.58 


25.32 


27.97 


30.53 


32.97 


-C{9',9') 


139.37 


62.18 


153.54 


22.58 


25.32 


27.97 


30.53 


32.97 


-0(10,10) 


139.37 


62.18 


153.54 


22.58 


25.32 


27.97 


30.53 


32.97 


-D(10',10') 


139.37 


62.18 


153.54 


22.58 


25.32 


27.97 


30.53 


32.97 


-0(7,7') 


-597.83 


-270.22 


-557.66 


-118.3 


-116.2 


-118.05 


-121. 


-124.28 


-0{7,9) 


641.09 


294.08 


582.32 


122.37 


126.07 


128.89 


131.1 


132.67 


-0(7,9') 


-232.82 


-105.66 


-228.83 


-42.16 


-44.54 


-46.71 


-48.75 


-50.66 


-0(7', 9) 


-246.68 


-112.02 


-241.42 


-44.81 


-47.26 


-49.49 


-51.58 


-53.53 


-0(7', 9') 


646.38 


296.49 


587.5 


123.33 


127.08 


129.95 


132.2 


133.82 


^(9.9') 


-110.27 


-50.13 


-107.04 


-20.14 


-21.19 


-22.13 


-23.01 


-23.81 


-0(10.1(1' ) 


-110.27 


-.^iD.i:-! 


-107.01 


-20.11 


-21.10 


-22,1:-! 


-2:i,01 


-23.81 



Table 9: Coefficients for the New Physics formula of {Fl). 
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(Pi) 




[1,6] 


[2,4.3] 


[4.3,8.68] 


[1,2] 


[2,3] 


[3,4] 


[4,5] 


[5,6] 


5 


0.051 


0.05 


0.059 


0.052 


0.056 


0.046 


0.046 


0.053 


^(0 0) 


-22.65 


-7.24 


-76.91 


0.85 


-1.45 


-4.02 


-7.1 


-10.93 




-62.37 


-25.4 


-87.61 


-12.98 


-10.73 


-11.09 


-12.59 


-14.98 


^(0 T) 


-781.92 


-286.75 


841.05 


-666.22 


-238.23 


-54.82 


51.6 


125.75 


^(0 9) 


-4.14 


-1.14 


-17.14 


0.31 


-0.13 


-0.68 


-1.38 


-2.27 


N(r) Qf\ 


-15.83 


-15.3 


124.84 


-19.74 


-12.07 


-3.79 


5.06 


14.72 


^(0 lO'l 


11.04 


4.2 


23.17 


0.73 


1.34 


2.07 


2.94 


3.96 




-136.39 


-56.11 


-223.52 


-11.23 


-19.14 


-27.21 


-35.35 


-43.45 


^(7 7) 


-3.85 


-3.33 


-77.22 


21.17 


3.71 


-4.55 


-9.97 


-14.21 


V ' : ' ) 


-312.25 


-136.11 


-200.72 


-95.67 


-66.13 


-54.77 


-49.26 


-46.42 


^(9 9) 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-^{9', 9') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-'^{10, 10) 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-^(10', 10') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 




-157.98 


-69.69 


-138.92 


-37.23 


-31.2 


-29.65 


-29.6 


-30.31 


^(7,9) 


425.78 


195.87 


378.22 


82.56 


84.39 


85.62 


86.41 


86.79 


-^^(7,9') 


-17.5 


-7.45 


-24.73 


-2.36 


-2.87 


-3.43 


-4.07 


-4.78 


^(7', 9) 


-26.64 


-11.65 


-32.85 


-4.13 


-4.68 


-5.27 


-5.92 


-6.64 


N(7' 9') 


426.35 


196.12 


378.93 


82.65 


84.49 


85.73 


86.54 


86.93 




-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-^{10,10') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


^(0 0) 


415.25 


143.29 


657.51 


108.6 


60.53 


61.39 


78.99 


105.74 


-^(0,7) 


-782.9 


-287.16 


8.38.78 


-666.19 


-238.35 


-55.03 


51.3 


125.37 


^(0 7') 


-45.61 


-19.25 


-105.64 


1.3 


-5.62 


-9.92 


-13.69 


-17.67 




-15.83 


-15.3 


124.84 


-19.74 


-12.07 


-3.79 


5.06 


14.72 


^(0 9') 


-4.14 


-1.14 


-17.14 


0.31 


-0.13 


-0.68 


-1.38 


-2.27 


^(n in"* 


-136.39 


-56.11 


-223.52 


-11.23 


-19.14 


-27.21 


-35.35 


-43.45 




11.04 


4.2 


23.17 


0.73 


1.34 


2.07 


2.94 


3.96 


^(7 7) 


7185.13 


3093.31 


2874.99 


2722.7 


1627.25 


1170.09 


915.02 


750.07 


\' ' ' J 


7195.21 


3097.72 


2882.27 


2725.55 


1629.34 


1171.9 


916.71 


751.71 


^(9,9) 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


D(9',9') 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 




31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


-0(10'. 10') 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


-D{7,7') 


-157.98 


-69.69 


-1.38.92 


-37.23 


-31.2 


-29.65 


-29.6 


-30.31 


■D(7,9) 


425.78 


195.87 


378.22 


82.56 


84.39 


85.62 


86.41 


86.79 


-D{7,9') 


-17.5 


-7.45 


-24.73 


-2.36 


-2.87 


-3.43 


-4.07 


-4.78 


-D{7',9) 


-26.64 


-11.65 


-32.85 


-4.13 


-4.68 


-5.27 


-5.92 


-6.64 


-D(7',9') 


426.35 


196.12 


378.93 


82.65 


84.49 


85.73 


86.54 


86.93 


-D(9.9') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-^''llO.Kj' ) 


-2.',(i 


-0.97 




-0.17 




-O.iS 


-O.CiS 


-0.92 



Table 10: Coefficients for the New Physics formula of (Pi). 
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{P2) 






To /I Q 1 


\ 4 Q Q fiQ 1 
[ 4.0 , O.DO J 






[ /1 1 


\ A c: 1 


[5,6J 


s 


0.066 


0.069 


0.064 


0.022 


0.031 


0.083 


0.082 


0.071 


^(0,0) 


34.88 


33.22 


-266.57 


42.57 


26.07 


8.3 


-10.68 


-31.39 


^(0,7) 


-917.02 


-421.92 


-813.51 


-177.91 


-181.82 


-184.41 


-186.07 


-186.81 


■'^(0,7') 


-13.39 


-6.16 


-11.88 


-2.6 


-2.66 


-2.69 


-2.72 


-2.73 


^(0,9) 


-67.97 


-27.97 


-111.15 


-5.6 


-9.55 


-13.57 


-17.62 


-21.63 


-'^(0,10) 


-8.1 


-7.71 


61.87 


-9.88 


-6.05 


-1.93 


2.48 


7.28 


-^(0,10') 


2.39 


1.11 


1.81 


0.48 


0.48 


0.48 


0.48 


0.46 


N(7,7') 


-157.98 


-69.69 


-138.92 


-37.23 


-31.2 


-29.65 


-29.6 


-30.31 


^(7,9) 


425.78 


195.87 


378.22 


82.56 


84.39 


85.62 


86.41 


86.79 


^{7,9') 


-17.5 


-7.45 


-24.73 


-2.36 


-2.87 


-3.43 


-4.07 


-4.78 


^(7', 9) 


-26.64 


-11.65 


-32.85 


-4.13 


-4.68 


-5.27 


-5.92 


-6.64 


A'(7',9') 


426.35 


196.12 


378.93 


82.65 


84.49 


85.73 


86.54 


86.93 


^{9,9') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-^(10,10') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-D{0,0) 


415.25 


143.29 


657.51 


108.6 


60.53 


61.39 


78.99 


105.74 


^(0,7) 


-782.9 


-287.16 


838.78 


-666.19 


-238.35 


-55.03 


51.3 


125.37 


-D(0,7') 


-45.61 


-19.25 


-105.64 


1.3 


-5.62 


-9.92 


-13.69 


-17.67 


-D{0,9) 


-15.83 


-15.3 


124.84 


-19.74 


-12.07 


-3.79 


5.06 


14.72 


■D(0,9') 


-4.14 


-1.14 


-17.14 


0.31 


-0.13 


-0.68 


-1.38 


-2.27 


-D(0,10) 


-136.39 


-56.11 


-223.52 


-11.23 


-19.14 


-27.21 


-35.35 


-43.45 


-'^{oao') 


1 1 r\A 




ZO.l ( 


n 

U. i o 


l.o4 


AMI 




o.yo 


-D(7,7) 


7185.13 


3093.31 


2874.99 


2722.7 


1627.25 


1170.09 


915.02 


750.07 


-D(7',7') 


7195.21 


3097.72 


2882.27 


2725.55 


1629.34 


1171.9 


916.71 


751.71 


-D{9,9) 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


^(9', 9') 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


-D(10,10) 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 




31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


-D(7,7') 


-157.98 


-69.69 


-138.92 


-37.23 


-31.2 


-29.65 


-29.6 


-30.31 


-C{7,9) 


425.78 


195.87 


378.22 


82.56 


84.39 


85.62 


86.41 


86.79 


■D(7,9') 


-17.5 


-7.45 


-24.73 


-2.36 


-2.87 


-3.43 


-4.07 


-4.78 


■D(7',9) 


-26.64 


-11.65 


-32.85 


-4.13 


-4.68 


-5.27 


-5.92 


-6.64 


-D(7',9') 


426.35 


196.12 


378.93 


82.65 


84.49 


85.73 


86.54 


86.93 


-D(9,9') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-D(10,10') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 



Table 11: Coefficients for tlie New Pfiysics formula of {P2)- 
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(Ps) 




[1,6] 


[2,4.3] 


[4.3,8.68] 


[1,2] 


[2,3] 


[3,4] 


[4,5] 


[5,6] 


S 


0.024 


0.024 


0.027 


0.026 


0.027 


0.022 


0.021 


0.024 




-1.43 


-0.6 


-0.48 


-0.4 


-0.28 


-0.24 


-0.24 


-0.27 


^(0,7) 


3.06 


1.27 


0.83 


1.06 


0.65 


0.49 


0.42 


0.43 


^(0,7') 


142.66 


59.32 


38.58 


49.47 


30.27 


22.92 


19.81 


20.19 


■'^(0,9) 


-0.13 


-0.06 


-0.06 


-0.02 


-0.02 


-0.02 


-0.03 


-0.03 


-^(0,9') 


8.97 


3.78 


3.99 


1.5 


1.57 


1.68 


1.87 


2.34 


■^(7,7') 


-157.98 


-69.69 


-138.92 


-37.23 


-31.2 


-29.65 


-29.6 


-30.31 


-^(7,9) 


425.78 


195.87 


378.22 


82.56 


84.39 


85.62 


86.41 


86.79 


^(7,9') 


-17.5 


-7.45 


-24.73 


-2.36 


-2.87 


-3.43 


-4.07 


-4.78 




-26.64 


-11.65 


-32.85 


-4.13 


-4.68 


-5.27 


-5.92 


-6.64 


^(7', 9') 


426.35 


196.12 


378.93 


82.65 


84.49 


85.73 


86.54 


86.93 


-^^(9,9') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-'^{10,10') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


^fO 01 


415.25 


143.29 


657.51 


108.6 


60.53 


61.39 


78.99 


105.74 


-^fO 71 


-782.9 


-287.16 


838.78 


-666.19 


-238.35 


-55.03 


51.3 


125.37 


P^(0 7') 


-45.61 


-19.25 


-105.64 


1.3 


-5.62 


-9.92 


-13.69 


-17.67 


^fO 91 


-15.83 


-15.3 


124.84 


-19.74 


-12.07 


-3.79 


5.06 


14.72 


^(0 9') 


-4.14 


-1.14 


-17.14 


0.31 


-0.13 


-0.68 


-1.38 


-2.27 


-D{o,io) 


-136.39 


-56.11 


-223.52 


-11.23 


-19.14 


-27.21 


-35.35 


-43.45 


■D(0,10') 


11.04 


4.2 


23.17 


0.73 


1.34 


2.07 


2.94 


3.96 


-D(7,7) 


7185.13 


3093.31 


2874.99 


2722.7 


1627.25 


1170.09 


915.02 


750.07 


-C(7',7') 


7195.21 


3097.72 


2882.27 


2725.55 


1629.34 


1171.9 


916.71 


751.71 


-D(9,9) 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


-D(9',9') 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


-D(10,10) 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


-'5(10', 10') 


31.66 


13.02 


51.88 


2.61 


4.44 


6.32 


8.21 


10.08 


-D{7,7') 


-157.98 


-69.69 


-138.92 


-37.23 


-31.2 


-29.65 


-29.6 


-30.31 


-D(7,9) 


425.78 


195.87 


378.22 


82.56 


84.39 


85.62 


86.41 


86.79 


-'5(7,9') 


-17.5 


-7.45 


-24.73 


-2.36 


-2.87 


-3.43 


-4.07 


-4.78 


-'5(7', 9) 


-26.64 


-11.65 


-32.85 


-4.13 


-4.68 


-5.27 


-5.92 


-6.64 


-^5(7', 9') 


426.35 


196.12 


378.93 


82.65 


84.49 


85.73 


86.54 


86.93 


-'5(9,9') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 


-'5(10,10') 


-2.56 


-0.97 


-5.38 


-0.17 


-0.31 


-0.48 


-0.68 


-0.92 



Table 12: Coefficients for the New Physics formula of (P3). 
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them in quadrature. In the case of the observables (Afb), (-^l), (-^1,2,3) and (P4 5g), the 
numbers C?^'^ and 5k are given by Eq. fl34|) . We add to the similar quadratic terms for 
the theoretical quantities involved (decay constants, form factors, quark masses, CKM 
matrix elements). 

We want to obtain the constraints in the two-dimensional plane corresponding to 
two (shifts of) Wilson coefficients ((5Ca, 5Cb) among all the parameters p (we denote the 
remaining theoretical parameters q). We obtain these regions by drawing the region of 
{SCa-i 5Cb) where: 

min x^(5Ca, 5Cb) — min < 5 (36) 

q p 

where 6 depends on the dimension of the parameter space where the region is drawn (here, 
a two-dimensional plane), and the required confidence level (here, 2.3 for 68.3% CL and 
5.99 for 95.5% CL, see the review of statistics and Monte Carlo techniques in Ref. |13] ) 

At the practical level, one could in principle compute the on a grid for the Wilson 
coefficients, performing the minimisation over all the nuisance parameters. One would 
obtain the two-dimensional contours for a given pair of Wilson coefficients {6Ca, SCb) by 
performing a further minimisation on the other Wilson coefficients. It turns out that 
one can improve the accuracy of the method by sampling the parameter space through 
a Metropolis-Hastings Markov-Chain Monte Carlo algorithm with a weight that favours 
the minima of the x^? for instance exp(— x^/2). Once this sampling has been performed, 
one creates the grid of points (and compute the value of the points) by considering all 
the points sampled by the algorithm. The points of grid which have not been sampled at 
all are given a very large x^- One then proceeds to a smearing procedure to get a smooth 
dependence of the reconstructed grid. 
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